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ABSTRACT 


One  of  the  primary  research  tools  in  high  energy  physics  is  the  collision  of  beams  of 
particles  with  stationary  targets.  The  useful  information  that  may  be  extracted  from  such 
experiments  is  proportional  to  knowledge  of  the  thickness  of  the  target.  This  determines 
how  many  point  constituents  w'ere  available  for  the  particles  in  the  beam  to  interact  with. 
It  is  possible  to  determine  the  thickness  of  a  target  with  statistical  errors  of  the  order  of 
one  percent  by  the  use  of  x-ray  attenuation  measurements. 

Six  target  samples  from  the  Stanford  Linear  Accelerator  Center  Experiment  #143, 
"Nucleon  Spin  Structure  at  30  GeV,"  were  measured  using  this  method.  The  targets  were 
frozen  ammonia  crystals.  A  40  mCurie  americium  241  source  was  used  to  pass  a 
collimated  beam  of  known  intensity  x-rays  through  the  samples.  The  x-ray  intensity  on 
the  other  side  was  measured  by  a  silicon  cell  detector.  The  ammonia  material  was 
removed  and  the  process  repeated.  By  comparing  the  intensities  with  and  without  the 
material  to  the  tabulated  attenuation  coefficients,  the  thickness  of  material  is  calculated. 
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I.  INTRODUCTION 


Nuclear  and  particle  physics  has  long  been  a  branch  of  science  that,  by  necessity,  takes  an 
indirect  approach.  Investigators  of  the  subatomic  world  developed  the  approach  of  accelerating  beams  of 
individual  charged  particles  at  bulk  collections  of  other  particles  (targets)  to  check  predictions  of 
theoretical  models.  In  191 1  Ernest  Rutherford  postulated  a  dense  nucleus  at  the  center  of  the  atom  in 
contrast  to  the  dispersed  “Rice  Pudding”  model  of  the  time,  because  his  scattering  experiments  sometimes 
gave  large  angle  scattering  events  that  could  only  be  explained  if  the  positive  charge  was  distributed  over 
a  relatively  much  smaller  volume  than  the  radii  of  the  “electron  cloud. 

To  probe  smallar  distance  scales  requires  beams  that  transfer  larger  momentum.  This  is  a 
consequence  of  the  Heisenberg  Uncertainty  Principle  which  states  that  the  product  of  the  root  mean  square 
uncertainties  in  position  and  momentum  is,  at  best,  of  the  order  of  Planck’s  Constant.  Therefore,  to  probe 
distances  ~lfin  requires  momentum  transfers  of  at  least  a  few  hundred  MeV/c.  By  the  mid-1960’s,  the 
Stanford  Linear  Accelerator  Center  in  Palo  Alto,  California  was  using  beams  of  electrons  of  such  high 
momentum  and  energy  in  their  experiments  that  they  began  to  observe  the  same  large  angle  scattering 
from  free  point-like  soiurces  within  individual  nucleons.  These  experimental  point-like  objects  were 
termed  “partons”  by  RP.  Feynman,  and  later  it  became  accepted  that  they  included  the  “quarks”  of  M. 
Gell-Mann.  Over  the  span  of  a  few  years,  the  theory  of  hadron  interactions  evolved  from  the  exchange  of 
Yukawa  mesons  between  iso-spin  singlets  of  SU(2)  (nucleons),  to  the  exchange  of  gluons  between  color 
singlets  of  SU(3)0  SU(3)  (quarks),  and  quantum  chromodynamics  was  bom. 

Over  two  decades  later,  the  theory  was  faced  with  what  became  known  as  “The  Proton  Spin 
Crisis.”  A  QCD  calculation  known  as  “The  Ellis-Jaffe  Sum  Rule”  predicted  that  the  spin  of  the  nucleon 
was  simply  the  sum  of  the  spins  of  the  constituent  quarks  (Ellis,  J.  and  Jaffe,  R,  1974).  Early  attempts  to 
verily  the  prediction  were  inconclusive  due  to  relatively  large  random  uncertainties  (Alguard,  M.J.,  et  ai, 
SLAC  E80, 1976  and  1978  and  Baum  G.  et  al,  SLAC  E130, 1983).  In  1988,  after  more  precise 
measurement  techniques  became  available,  it  was  reported  by  experimenters  at  CERN  that  their 
observation  of  scattering  events  did  not  support  the  Ellis-Jaffe  prediction  (Ashman,  J.,etal.,  EMC,  1988 
and  1989).  Over  the  next  six  years,  collaborations  of  physicists  hurried  to  repeat  this  experiment  and 
check  the  findings  both  at  CERN  and  SLAC  (Abe,  K.,  et  al.,  1994,  Adams,  D.,  et  al.,  SMC,  1994,  Adeva, 
B.,  et  al..  SMC,  1993  and  Anthony,  P.L.,  et  al.,  E142,  1993). 

The  most  recent  of  these  was  the  E143  experiment.  E143  used  polarized  electrons  of  up  to  29. 1 
GeV  to  bombard  a  polarized  ammonia  target.  In  order  to  determine  the  number  of  protons  available  for 
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the  beam  to  interact  with,  and  thus  evaluate  the  desired  experimental  quantities,  it  was  necessarj'  to 
determine  the  thickness  of  the  target.  This  thesis  is  the  report  of  the  measurement  of  the  E143  target 
inserts  by  comparing  relative  intensity  data  from  the  attenuation  of  Am  x-rays  passing  through  the 
target  material.  The  advantage  of  the  E143  setup  was  that  using  electrons  at  SLAC  (as  opposes  to  the 
higher  energy  muons  at  CERN)  allowed  for  random  beam  polarization  flips  and  much  higher  currents 
and,  consequently,  event  rates.  These  two  factors  allowed  for  far  better  experimental  statistics  in  this 
energy  range  than  ever  aceived  before  (Abe,  K.,  et  ai,  1994  and  McCarthy,  J.,  et  al,  1993).  In  order  to 
preserve  the  advantage  of  low  experimental  error  in  the  experiment  overall,  the  contributions  of  error 
from  other  variables  such  as  beam  polarization,  target  polarization  and  target  thickness  had  to  be 
minimized.  The  value  of  the  target  thickness  needed  to  be  known  to  better  than  3%  fractional  uncertainty 
in  order  to  preserve  the  desired  low  error  of  the  entire  experiment  (McCarthy,  J.,  et  ai,  1993).  This  thesis 
IS  a  record  of  the  experimental  work  that  was  performed  by  myself  and  my  advisor  from  September,  1993 
through  September,  1994  to  precisely  and  accurately  measure  the  ammonia  target  thicknesses  used  in  the 
Stanford  Linear  Accelerator  Center  Experiment  #143.  Reduction  of  the  random  uncertainty  of  this 
experiment  directly  contributes  to  international  understanding  of  what  has  come  to  be  known  as  the 
“Nucleon  Spin  Crisis.” 
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II.  GENERAL  THEORETICAL  BACKGROUND 


The  modem  theory  of  the  strong  interaction  between  hadrons  is  called  quantum  chromodynamics  or 
“QCD.”  As  its  nama  implies  the  prefix  “chromo-”  refers  to  color.  In  the  relativistic  quantum  field  theory 
required  to  describe  QCD,  the  quantum  “numbers”  assigned  to  the  individual  constituents  are  called 
“color.”  The  individual  constituents  of  QCD  theory  are  called  “quarks”  and  they  are  held  together  in  the 
form  of  hadrons  via  the  color  fields  mediated  by  “gluons.”  This  is  analogous  to  the  individual  constituents 
of  QED  (quantum  electrodynamics),  the  electrons,  being  held  together  in  patterns  around  the  nucleus  via 
the  Maxwell  fields  mediated  by  photons.  There  are,  however,  several  critical  differences  between  the  well 
known  QED  and  the  newer  QCD  theories.  The  color  fields  of  QCD  theory  are  examples  of  non-Abelian 
gauge  fields  while  the  Maxwell  fields  of  QED  are  Abelian  gauge  fields.  The  difference  is  that  the  force 
mediators,  bosons,  of  an  Abelian  gauge  theory  are  not  influenced  by  the  force  they  carry ,  i.e.  they  do  not 
carry  the  “charge”  of  their  force.  In  a  non-Abelian  gauge  field  theory,  however,  the  force  mediating 
bosons  do  feel  the  strong  force  and  do  carry  the  color  charge  themselves.  This  leads  to  very  powerful,  non 
linear  self  interactions  and  so  the  conventional  perturbation  theory  of  canonical  quantum  mechamcs  is 
invalid.  Two  more  remarkable  results  of  this  decidedly  different  mathematical  description  of  the  strong 
force  are  “color  confinement”  and  “asymptotic  freedom.”  Color  confinement  refers  to  the  fact  that  all 
physical  bound  hadron  states  are  composed  of  multiple  quarks  whose  net  color  charge  sums  to  zero,  i.e., 
lone  quarks  with  a  “bare”  color  charge  are  not  observed.  Asymptotic  freedom  means  that  the  closer 
together  quarks  are,  the  less  they  feel  the  strong  force.  In  the  limit  of  zero  range  (requiring  infinite  energy 
by  the  Uncertainty  Principle  mentioned  above),  they  are  “free”  particles  and  feel  no  strong  force  at  all. 

One  way  to  explain  these  two  phenomena  is  to  consider  the  partons/quarks  as  eigenstates  of  the  free  field 
Hamiltonian  and  thus  only  truly  realizable  in  the  limit  of  infinite  energy  (Feynman,  R.P.,  1972). 

The  most  enduring  model  of  hadron  structure  is  the  “naive”  quark-parton  model  (Bjorken,  J.D.,  1966 
and  Feymnan,  R.P.,  1972).  In  this  model,  the  nucleon  is  composed  of  three  “constituent”  quarks  of 
various  “flavors,”  primarily  “up”  and  “down”  for  the  vast  majority  of  visible  matter  in  the  universe.  The 
flavors  are  quantum  numbers  of  the  global  SU(3)  symmetry  of  the  strong  force,  while  the  colors  are 
quantum  numbers  of  the  local  SU(3)  symmetry.  The  Uncertainty  Principle  accounts  for  yet  another 
strange  effect.  On  the  sub  fermi  scale,  the  “vacuum”  in  relativistic  quantum  field  theory  is  constantly 
creating  and  annihilating  particles  and  anti-particles  with  masses  (energies)  inversely  proportional  to  their 
lifetimes.  This  means  that  the  constituent  quarks  inside  of  hadrons  are  also  immersed  in  a  “virtual  sea”  of 
other  particles,  primarily  “strange”  quarks.  This  constituent  quark  model  has  been  very  successful  in 
describing  global  hadron  properties  such  as  charge  and  mass  so  it  was  assumed  that  it  would  work  equally 
well  in  predicting  the  spin  of  the  hadron  to  be  the  superposition  of  the  spins  of  the  effective  fields  of  the 
constituent  quarks.  Surprisingly,  this  was  not  the  case  (Ashman  et  al.,  EMC,  1988  and  1989). 
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A.  DEEP  INELASTIC  E-P  SCATTERING 


One  of  the  primaiy  modem  tools  of  testing  QCD  predictions  is  deep  inelastic  scattering  of  leptons 
by  hadrons.  There  are  several  reasons  that  this  is  such  a  desirable  method  of  probing  nucleon  structure. 
First  of  all,  leptons  do  not  experience  the  strong  force  and  so  the  reaction  is  governed,  to  first  order,  by 
single  interaction  kinematics.  Secondly,  since  leptons  have  such  small  masses,  they  may  be  accelerated  to 
an  ultra  relativistic  state  by  lesser  applied  forces  than  those  required  to  produce  an  equal  effect  in  more 
massive  probes.  The  calculations  are  performed  in  the  limit  of  infinite  momentum.  What  this  means  is 
that  the  state  of  asymptotic  freedom  mentioned  above  is  assumed  to  be  in  effect.  Tests  of  the 
electromagnetic  stmcture  fimctions  of  the  nucleon  provide  insight  into  the  parton  distribution  inside.  The 
name  deep  inelastic,”  means  that  a  great  deal  of  the  incident  electron’s  momentum  is  transferred  to  the 
target  nucleon,  and  so  a  finer  structure  may  be  discerned  as  discussed  in  the  introduction.  Since  modem 
electron  beams  are  in  the  ultra  relativistic  regime  (y  ~  lO')  and  nucleon  targets  are,  by  definition,  at  rest  in 
the  lab  frame  of  reference,  it  is  convenient  to  analyze  the  kinematics  in  that  frame.  Figure  (1)  provides 
an  illustration  of  the  relation  between  the  beam  axis,  taken  to  be  in  the  +z  direction,  and  the  nucleon, 
taken  to  be  at  rest  in  the  laboratory  reference  frame. 


p  -  cp**-  (E,  cp)  =  electron  four-momentum^/Rnal  electron  state:|p') 
Initial  electron  state:  |p,  \  Scattering  angle:  0 


Virtual  photon  momentum:  q,  e  y 


linalji^Ieon  slates: 


Initial  nucleon  state:  P  -  cP"  =  (M,  0)  =  nucleon  four-momentum 

9^  =  q\  =  -^  =  (pV’^=  -4EE’Sin^(e/2)  =  -2EE’[1-Cos(0)]  <  0 


Figure  1.  Schematic  diagram  of  polarized  deep  inelastic  electron-nucleon  scattering. 

Other  kinematic  variables  used  to  describe  the  collision  interaction  are  (Bjorken,  J.D.,  1970  and  Bjorken, 
J.D.  and  Paschos,  E.A.,  1969,  Feynman,  R.P.,  p.  126,  1972,  Hughes,  V.W.  and  Kuti,  J.,  1983  and  Kaku 
M.,  1993,  pp.  460-461): 

E  =  Energy,  in  GeV,  of  the  incident  electron  beam  in  the  lab  reference  frame. 

E’  =  Energy,  in  GeV.  of  the  scattered  electron  beam  in  the  lab  reference  frame 
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m  =  nieC^  =  rest  mass  energy  of  the  electron  ~  0.51  MeV. 

M  =  MpC^  =  rest  mass  energy  of  the  nucleon  ~  940  MeV  =  (P^QJ/(E-E’). 

V  =  E-E’  =  (FQj/A/=  energy  transferred  from  electron  to  nucleon  by  virtual  photon. 

E  =  {l  +  2j^l  +  Tan^  (%)}  =  polarization  of  virtual  photon. 

s  =  »»'*(cp,  0,  0,  +/-  E)  =  spin  polarization  four-vector  of  the  electron. 

S  =  (0,  o)  =  spin  polarization  four-vector  of  the  nucleon. 

=>  =  p*‘p^  =  5^  =  S-S^  =  -1,  P^S,  =  0. 

B.  ELECTROMAGNETIC  FORM  FACTORS 

Even  though  lepton-hadron  scattering  is  a  strictly  electromagnetic  process  at  accelerator  energies,  the 
strong  interaction  modifies  the  hadron  electromagnetic  current  at  the  sub-fermi  size  scale.  In  the  “naive” 
quark/parton  model,  the  interaction  is  assumed  to  occur  due  to  hard  scattering  of  a  virtual  photon  off  of  a 
single  constituent  of  the  nucleon,  possessing  some  positive,  real  fraction  of  the  total  nucleon  spin  and 
momentum.  The  differential  cross  section  for  deep  inelastic  scattering  of  polarized  electrons  from 
polarized  nucleons  in  the  first  Bom  approximation,  with  factors  of  Planck’s  constant  suppressed,  is 
(Hughes,  V.W.  and  Kuti,  J.,  1983,  Kaku,  M.,  1993,  p.  462  and  McCarthy,  J.,  et  al.,  E143,  1993) 

d^dE'  le'J  IeJ  ''''  4EE'  +  (2' 

In  this  equation,  is  the  charged  leptonic  tensor  amplitude,  H'”  is  the  charged  hadronic  tensor 
amplitude,  a  is  the  electromagnetic  fine  stmctme  constant,  (E’/E)  is  the  relativistic  recoil  effect  term  and 
Ojwon  is  the  Mott  scattering  cross  section  (Bosted,  P.,  1993), 

■2aE'Cos(^)T  r  aCos{%)  f 
q2  j  -[2ESin^(?^)J  ■ 

The  leptonic  tensor  amplitude  is  defined  as  the  contribution  of  the  electron-photon  vertex  in  the  Feynman 
diagram  for  the  scattering  process.  It  may  be  decomposed  into  a  spin-independent  symmetric  part  and  a 
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spin-dependent  anti-symmetric  portion  (Feynman,  R.P.,  p.  126,  1972,  Hughes,  V.W.  and  Kuti,  J.,  1983 
and  Kaku,  M.,  1993,  p.  462).  Summed  over  final  polarization  states  of  the  scattered  electron,  becomes 

=  (p'p^  -1-  p;p^  -g;,vP"Pa)  +  ^("ie^vapS“q^),  (3] 

where  s!"'^  is  the  completely  anti-symmetric  fourth  rank  tensor,  e'”''=l,  and  is  the  Minkowski  metric. 
In  a  similar  fashion,  the  hadron  tensor  may  also  be  decomposed  into  a  spin-independent  symmetric  part 
and  a  spin-dependent  anti-symmetric  portion  (Ellis,  J.  and  Jaffe,  R.,  1974,  Feynman,  R.P.,  p.  126,  1972, 
Hey,  A.J.G.  and  Mandula,  J.E.,  1972,  Hughes,  V.W.  and  Kuti,  J.,  1983,  Kaku,  M.,  1993,  pp.  462,  Kuti,  J. 
and  Weisskopf,  F.,  1971  and  Pokorski,  1987,  p.  176).  Written  in  terms  of  the  expectation  value  of  current 
operators  between  initial  and  final  states  H*’  becomes 

i 

=  '^{27ty5{q  +  P-P'){UP,  5|/(o,o)f/(o)|/){/|f/-^(0);''(o,o)|f/P,  S),  (4) 

where  J**  is  the  fully  interacting,  Heisenberg  picture,  hadronic  current  operator  for  lepton  scattering  in  the 
global  SU(3)  approximation  and/  is  the  corresponding  Schroedinger  picture  “free”  or  “bare”  current.  It 
can  be  shown  (Drell,  S.D.,  Levy,  D.J.  and  Van,  T.,  1969  and  1970,  Feynman,  R.P.,  1972,  Ji,  X.,  1993, 
Kaku,  M.,  1993  and  Kuti,  J.  and  Weisskopf,  F.,  1971)  that 7“  is  also  the  sum  of  a  symmetric,  conserved, 
vector  part  and  an  anti-symmetric,  partially  conserved  axial  part.  The  sum  over  ’  in  Equation  (4)  is 
over  all  possible  final  hadron  states.  After  the  smnmation,  the  expression  for  the  hadronic  tensor 
amplitude  further  reduces  to  (Drell,  S.D.,  Levy,  D.J.  and  Van,  T.,  1969  and  1970,  Hughes,  V.W.  and  Kuti, 
J.,  1983,  Kaku,  M.,  1993,  p.  462  and  Pokorski,  1987,  p.  175) 

=  5|[y'(O.J).J'(o.o)]|/>.  5).  ,5, 

In  this  step,  the  operator  product  of  the  currents  is  converted  to  a  commutator  because  the  two  differ  by  a 
term  proportional  to  6(q  +  p  -  p’)  which  is  always  identically  zero  in  physical  situaUons  where  q  and  Ap 
are  always  strictly  positive  (Kaku,  M.,  1993,  pp.  462  and  Pokorski,  1987,  pp.  175).  The  commutator  may 
then  be  evaluated  using  light  cone  analysis  Prell,  S.D.,  Levy,  D.J.  and  Van,  T.,  1969  and  1970  and 
Feynman,  1972)  to  give  a  theoretical  value  to  compare  to. 
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The  spin-independent  symmetric  part  of  the  hadron  tensor  may  be  completely  defined  by  two  real 
scalar  spin-independent  structure  functions,  IF]  and  Wj.  From  conserved  parity  and  current  arguments,  as 
well  as  Lorentz  and  gauge  invariance,  the  symmetric  part  of  the  hadronic  tensor  amplitude  must  have  the 
general  form  (Drell,  S.D.,  Levy,  D.J.  and  Yan,  T.,  1969  and  1970,  Feynman,  R.P.,  p.  127,  1972,  Hey, 
A.J.G.  and  Mandula,  J.E.,  1972,  Jaffe,  R.L.,  1975,  Ji,  X.,  1993  and  Kaku,  M.,  1993,  p.  462) 


TT  MV  _ 
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V  Q"  j 


(6) 


averaged  over  initial  polarization  states  of  the  nucleon.  Also  note  that  definitions  of  structure  fimctions 
throughout  the  references  differ  by  factors  of  A/^/v  used  for  normalization  (Bosted,  P.,  1993  and  Ji,  X., 
1993).  Similarly  an  expression  for  the  spin-dependent,  anti-symmetric  part  of  the  hadronic  tensor 
amplitude  may  be  given  in  terms  of  two  spin-dependent  structure  fimctions,  G,  and  G2  (Hey,  A.J.G.  and 
Mandula,  J.E.,  1972,  Hughes,  V.W.  and  Kuti,  J.,  1983  and  Ji,  X.,  1993)  and  H(Ar  becomes 

H„;"  =  £"«'q.{A/S,G,+L[(p^,‘)Sj-(5j,s)Pj]Gj,  (7) 


1.  The  Bjorken  scaling  limit 

Bjorken  theorized  (Bjorken,  J.D.,  1966  and  1970)  that  in  the  scaling  limit  as  v  and  <2^  go  to 
infinity,  the  structure  functions  become  dependent  not  on  v  and  but  rather  on  their  ratio  and  the  mass 
of  the  target  nucleon.  This  ratio  became  known  as  the  scaling  variable  “x-BJorken”  =  ^IlMv.  In  his 
original  papers  (Bjorken,  J.D.,  1966  and  1970),  he  used  the  symbol  to  which  is  defined  as  the  inverse  of 
XBj.  In  the  Bjorken  scaling  limit  (infinite  momentum  frame  =>  asymptotic  freedom  of  quarks),  where  v 
and  ^  go  to  infinity  with  Xbj  held  fixed,  both  the  spin  dependent  and  the  spin  independent  structure 
functions  of  v  and  scale  to  new  “asymptotic”  structure  functions  which  depend  only  on  xbj.  The 
references  vary  on  convention  of  definitions  in  terms  of  factors  ofA/^/v  (Bosted,  P.,  1993  and  Ji,  X., 
1993).  W\  and  W2  scale  to  F\  and  Fi  where  F,  is  the  “Dirac”  form  factor  which  reduces  to  the  nucleon 
electric  charge  in  the  limit  as  0  and  F2  is  the  “Pauli”  form  factor  which  reduces  to  the  nucleon 
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anomalous  magnetic  moment.  In  a  similar  fashion,  the  spin-dependent  structure  functions  G|  and  G2 
become  the  asymptotic  spin-dependent  structure  functions  gi  and  g2. 

The  deep  inelastic  form  factors  may  be  studied  using  the  techniques  of  current  algebra’s 
evaluated  on  the  light  cone  defined  by  x\  =  0.  The  theoretical  mathematical  techniques  of  evaluating 
current  commutators  in  the  infinite  momentum  approximation  were  pioneered  by  M.  Gell-Mann  of  the 
California  Institute  of  Technology  and  J.  Schwinger  of  Harvard  University.  One  valuable  experimental 
prediction  using  these  techniques  was  made  by  J.D.  Bjorken  at  SLAC.  He  predicted  the  behavior  of  the 
difference  of  the  first  asymptotic  structure  functions  of  the  neutron  and  proton.  The  Bjorken  sum  rule 
prediction,  with  third  order  QCD  corrections  (K.  Abe  et  al.,  1994,  B.  Adeva  et  ai,  SMC,  1993,  D. 
Adams  et  al.,  SMC,  1994,  Bjorken,  J.D.,  1966  and  1970  and  Feynman,  R.P.,  1972),  is 
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The  term  in  parenthesis  on  the  right  is  a  radiative  correction  where  Ok{^)  is  the  energy-dependent  strong 
coupling  constant  and  gA  and  gv  are  the  vector  and  axial-vector  weak  coupling  constants  from  nucleon 
beta  decay  data.  The  ratio  of  gA  to  gv  is  found  experimentally  to  be  1.257  +/-  0.003. 

Bjorken  used  strict  U(6)  local  and  global  symmetry  assumptions  to  derive  his  rule.  This  is  the  strictest 
treatment  of  quarks  as  it  considers  all  6  types  known  in  the  standard  model.  The  vast  majority  of  all 
physical  processes,  however,  only  involves  the  3  types  already  mentioned,  up,  down  and  strange,  and  so 
dealing  with  only  those  three  flavors  of  quarks  should  be  an  adequate  approximation.  This  approximation 
was  used  by  J.  Ellis  and  R.  Jaffe  to  make  a  prediction  of  the  separate  first  asymptotic  structure  functions  of 
the  neutron  and  proton.  The  Ellis-Jaffe  sum  rule  prediction,  under  the  assumptions  of  global  SU(3)  flavor 
symmetry  and  no  net  polarization  of  the  (strange)  sea  quarks,  predicts  (K.  Abe  et  al.,  1994,  B.  Adeva  et 
ai,  SMC,  1993,  D.  Adams  et  ai,  SMC,  1994,  Ellis,  J.  and  Jaffe,  R.,  1974)  a  value  for  g.P,  integrated  over 
values  of  Bjorken  x  frqm  zero  to  one  of  0. 160  +/-  0.006.  The  failure  of  the  experimental  evidence  to 
support  this  prediction  has  caused  great  volumes  of  theoretical  debate  concerning  the  validity  of  the 
assumptions  and  possible  explanations  for  the  experimental  results  being  lower  than  theory. 

2.  SLAC  Experiment  #143:  “Nucleon  Form  Factors  at  30  GeV” 


The  E143  “Nucleon  Spin  Structure”  experiment  used  a  polarized  solid  ammonia  target  as  a 
source  of  polarized  proton  and  deuteron  scattering  centers.  The  target  was  placed  in  the  beamline  at  End 
Station  'A'.  The  electrons  were  accelerated  to  energies  of  29.1  GeV.  The  electrons  were  polarized  at  the 
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source  by  a  laser  striking  a  strained  GaAs  epitaxial  layer  on  the  cathode  and  imparting  the  laser 
polarization  to  the  liberated  electrons  (Abe,  et  aL,  1994  and  McC^arthy,  J,,  et  al.,  1993).  The  differential 
cross  section  for  unpolarized  lepton-hadron  scattering,  in  the  approximation  of  zero  lepton  mass  {M»m\ 
is  (Feynman,  R.P.,  p.  128,  1972) 


dadE' 


yQj 


(9) 


For  the  large  scattering  angles  of  interest,  this  differential  cross  section  is  effectively  equal  to  the 
unpolarized  cross  section,  Ounp  (McCarthy,  J.,  et  aL,  p.  3,  1993).  The  sign  of  the  anti-symmetric  hadron 
tensor  in  Equation  (1)  depends  on  the  helicity  of  the  proton  in  the  center  of  mass  reference  frame.  This 
means  that  the  sum  or  difference  of  the  cross  sections  for  two  scattering  angles  rotated  by  n  radians  will 
cancel  either  the  symmetric  or  anti-symmetric  part,  respectively.  The  sum  of  the  cross  sections  is, 
therefore,  simply  equal  to  twice  Cunp.  The  difference  of  the  cross  sections,  however,  is  the  experimental 
quantity  of  interest  because  it  has  a  definite  sensitivity  to  the  initial  nucleon  polarization.  It  is  (Ashman, 
J.,  et  al,  p.  3,  1989,  Hey,  A.J.G.  and  Mandula,  J.E.,  p.2610,  1972,  Hughes,  V.W.  and  Kuti,  I,  p.  216, 
1983,  McCarthy,  J.,  et  aL,  p.  3,  1993) 
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The  longitudinal  asymmetry  of  scattering  events,  A\oog,  is  then  given  in  terms  of  these  sums  and 
differences  of  cross  sections  by  (Abe,  et  al.,  1994,  Bosted,  1994) 
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where  Pb  represents  the  beam  polarization,  Pj  represents  the  target  polarization,  Cn  is  a  correction  for  the 
effects  of  polarization  of  the  unpaired  proton  in  A  is  the  munber  of  scattered  electrons  per  incident 
charge  for  negative  (N.)  and  positive  (A+)  beam  helicity,  Aocp  is  the  obseiv’ed  experimental  asymmetry  ,  /Irc 
is  the  QED  radiative  correction  for  higher  order  terms  of  the  fine  structure  constant  and  Fd  is  the  ratio  of 
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scattering  events  originating  from  polarizable  protons  to  the  total  number  of  scattering  events  defined  as 
(Ashman,  J.,  et  al.,  p.  12,  1989  and  Rondon,  O.A.,  UVA-INPP-PUB-91-15,  1991)  A^pap/  SyA^yOy,  where  Ap 
is  the  number  of  polarizable  protons  per  square  centimeter,  Cp  is  the  cross  section  for  electron  scattering 
from  protons,  and  the  summation  in  the  denominator  is  over  all  possible  scattering  sources.  A  similar 
expression  may  be  derived  for  the  perpendicular  asymmetry  of  scattering  events,  ^perp.  The  asymptotic, 
spin-dependent  structure  functions  may  then  be  derived  from  expressions  in  terms  of  the  cross  sections, 
kinematic  variables  and  Bjorken  x  (Abe,  K.,  et  al.,  E143,  1994,  Bosted,  1993  and  McCarthy,  J.,  et  al.). 

The  target  used  in  End  Station  A’  (ESA)  at  SLAC  for  E143  contained  grains  of  frozen  ammonia  that 
had  been  “radiation  damaged”  to  produce  “paramagnetic  impurities”  at  various  electron  beam  facilities 
including  the  Massachusetts  Institute  of  Technology  (MIT)  Bates  Laboratory,  Stanford  High  Energy 
Physics  Lab  (HEPL),  Naval  Postgraduate  School  (NPS)  and  Saskatchewan  Accelerator  Laboratory  (SAL). 
The  ammonia  in  the  target  is  bombarded  with  microwaves  at  a  frequency  to  cause  spin-flip  transitions  to 
polarized  states  around  the  impurities.  The  target  material,  target  housing  assembly  and  target 
polarization  method  are  described  in  depth  by  the  target  group  in  various  references  (Abe,  K.,  et  al.,  E143, 
1994,  Crabb,  D.G.,  1990  and  McCIarthy,  J.,  et  al.).  The  basic  experimental  setup  of  E143  in  ESA  is 
illustrated  in  Figure  2  below. 


Figure  2.  E143  experiment  schematic.  Polarized  electrons  are  scattered  by  a  polarized  ammonia  target. 
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III.  TARGET  MATERIAL  THICKNESS  MEASUREMENT 


The  ammonia  target  thickness  measurements  were  made  at  SLAC  from  September  1993  through 
February  1994.  The  work  was  performed  by  the  E143  Collaboration  target  group  including  myself  and  my 
advisor.  The  thickness  was  computed  using  the  principle  of  attenuation  of  radiation  through  matter.  A 
source  of  known  intensity  emits  x-rays  that  pass  through  the  unknown  quantity  of  ammonia  and  the 
reduced  intensity  is  detected  on  the  other  side.  By  comparison  of  the  emitted  and  detected  intensities,  the 
amount  of  material  passed  through  may  be  determined  to  a  high  degree  of  precision.  The  measurement 
apparatus  was  located  in  the  small  building  adjacent  to  the  controlled  entrance  to  End  Station  A  so  that 
minimum  changeout  time  and  effort  was  required  to  transport  the  targets  to  and  from  the  beamline  and 
thickness  measurement  area.  The  torlon  cells  in  the  target  insert  sticks  containing  the  ammonia  material 
were  perforated  to  allow  the  liquid  helium  bath  to  permeate  and  adequately  cool  the  cell  and  yet  not  allow 
the  beads  of  material  that  were  jostled  around  in  transit  to  be  shaken  out  of  the  cell. 


A.  APPARATUS  CONFIGURATION 

A  schematic  representation  of  the  liquid  helium  refrigerant  dewar  that  contained  the  target  insert 
and  ammonia  material  as  well  as  the  locations  of  the  quantities  of  interest  is  given  in  Figure  3  below. 


frozen  ammonia  material 
of  specified  density  and 
attenuation  coefficient 


incident  x-rays 


> 


liquid  helium  bath 
surroimding  and 
permeating  cells 


ciyostat  apparatus  including 
target  cells,  NMR  coils,  RF 
wave  guides,  end  caps,  etc. 


attenuated  x-rays 


Figure  3.  Schematic  of  target  thickness  determination  method. 
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The  following  is  a  list  of  symbolic  notation  for  the  quantities  used  in  the  measurement: 


ty=  effective  linear  thickness  of  material  “y”  in  cm 
ftrgt=  target  cell  measured  value:  3.00  +/-  0.03  cm 
tiTom—  effective  linear  thickness  of  ammonia  in  target  cell 

py  =  mass  density  of  material  "y"  in 

cm 

Xy  =Pyty  =  mass  thickness  of  material  "y"  in  — 

cm 

=  thickness  of  '‘He  outside  target  cell 
jc'hc  =  thickness  of  '*He  inside  target  cell 

=  thickness  of  cryostat,  target  assembly,  etc. 

=  mass  thickness  of  ammonia  in  target  cell 


Py  =  linear  attenuation  coefficient  of  "y"  in  cm  ' 


pj-  -  attenuation  coefficient  of "  y"  in 
fy  =  volume  packing  fraction  of  material  "y" 

Xmm  =  ammonia  fraction  of  target  volume  =  = 


cm 

gram 


•tTgl 


Pamm^I 


trgt 


The  mass-attenuation  coefficients  are  calculated  by  the  program  “PHOTCOEF’  (AIC  software,  1993) 
running  on  a  personal  computer  with  a  386  processor.  The  entering  arguments  for  the  program  are  the 
chemical  composition  of  the  target  material  and  the  energy  of  the  incident  x-rays.  For  the  40  milli-Curie 
Americium  241  source  used,  the  x-ray  energy  was  59.53  keV.  This  gives  a  computed  value  for  the  mass 
attenuation  coefficient  for  ^He  of 


,  =0.16534  -22^ 

<He  gram 


(12) 


PHOTCOEF  gives  value  for  “normal”  chemical  compounds  only  and  a  correction  is  required  for 
ammonia  of  different  isotopic  compositions.  The  photon  interaction  is  an  electromagnetic  process 
governed  by  the  charge  distribution  of  a  scattering  center  (molecule  in  this  case).  The  actual  values  of  the 
mass-attenuation  coefficients,  however,  are  derived  in  terms  of  interaction  cross  section  per  unit  mass 
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cross  section 
unit  mass 


(13) 


where  A  is  the  molecular  weight,  n  =  Navp/A  is  the  number  of  scattering  centers  per  unit  volume  and  Nav 
is  Avogadro’s  number.  Expanding  upon  this  definition,  we  derive  an  expression  for  the  mass-attenuation 
coefficients  for  and  '^NDs  in  terms  of  the  mass-attenuation  coefficient  for  '''NH3. 


Note  that  using  this  derivation  and  the  definitions  at  the  beginning  of  this  Section,  it  is  observed  that 
the  expression  for  the  Dilution  Factor  used  in  the  expression  for  ^4  long  can  also  be  written  as: 


We  assume  that  the  x-ray  scattering  amplitudes  is  proportional  to  the  charge  distribution  of  the  scattering 
center,  independent  of  its  mass,  i.e.  Oisnhs  =  Oi4nh3.  With  this  assumption,  the  expression  for  (p/p)  for 
'^NHs  in  terms  of  (p/p)  for  '''NH3  from  Equation  (14)  can  be  simplified. 


Note  that  it  is  the  ratio  of  the  molecular  weights  that  determines  the  isotopic  mass-attenuation  coefficient. 
At  59.53  keV,  the  coefficients  for  “normal”  '''NH3  and  the  isotopes  used  in  the  experiment  are 


^  =  (%)(;^),«.  .  07) 

=0->«836i  =  {%){%), .^_ 
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where  the  isotopes  are  corrected  by  the  ratio  of  the  number  of  nucleons.  Binding  energy  correcUons  are 
negligible,  amounting  to  less  than  0.05%  (Walker,  F.W.,  1984).  The  mass  density  values  used  are 

Ph,  =  0.128  ±  0.003 

p„,  =  0.917  ±  0.001  .  (18) 

p„,  =  1.056  ±  0.001  .*1= 

cm  _ 

All  densities  are  for  4.0  +/-  0. 1  Kelvin  and  liquid  helium  vapor  pressure  (Anderson,  H.J.,  1989).  The 
uncertainty  in  the  temperature  is  a  conservative  estimate  and  the  densities  apply  within  the  range. 

There  are  two  experimental  configurations  which  must  be  considered.  In  the  first  setup,  the  cryostat 
assembly  is  full  of  helium  at  4  K  with  a  target  stick  inserted.  The  x-ray  beam  is  scanned  across  the  face  of 
a  cell  that  contains  the  unknown  thickness  of  solid  ammonia.  The  second  setup  is  identical  except  that  the 
ammonia  crystals  have  been  removed  and  the  target  cell  is  empty.  The  exposure  time  for  both 
configurations  is  the  same,  so  Poisson  statistics  apply  (Taylor,  J.R.,  1982).  The  resultant  x-ray  intensities, 
“I  m,”  and  “lempty,”  after  traversing  the  targets  are; 


(19) 

for  the  x-rays  passing  through  an  ammonia  filled  cell  and 


for  the  x-rays  passing  through  an  empty  cell.  The  terms  with  the  subscript  “cryo-”  refer  to  the  effects  of 
the  experimental  apparatus  and  include  all  support  structure,  NMR  coils,  target  cells,  etc.  These  are  too 

complex  to  be  easily  evaluated  but  they  are  identical  for  both  situations.  Consequently,  taking  the  ratio 
cancels  their  effects. 
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Taking  the  logarithm  of  each  side  and  rearranging  terms  gives  the  desired  expression  for  the  volume 
packing  fraction  in  a  target  cell. 


The  denominator  in  this  expression  can  be  computed  independently  of  the  x-ray  measurement  and  has 
values  of  0.47689  for  the  NH3  and  0.46988  for  the  ND3.  The  mass-thickness  of  ammonia  in  a  target  cell 
is  consequently: 


(23) 


The  denominator  has  a  value  of  0. 17335  for  the  NH3  and  0. 14832  for  the  ND3.  Note  that  the  effective 
mass-attenuation  coefficient  of  helium  in  the  denominator  is  suppressed  relative  to  the  mass-attenuation 
coefficient  of  ammonia  by  a  factor  of  the  ratio  of  their  densities  (<  0. 143).  Thus  while  each  mass- 
attenuation  coefficients  is  comparable  in  magnitude,  the  ammonia  term  contributes  (1/0.143)  =  7  times 
more  in  the  denominator.  This  observation  will  have  implications  in  the  next  section. 


B.  EXPERIMENTAL  PROCEDURES 


We  assume  that  if  x,...,  z  are  measurements  used  to  compute  the  function  q(x,...,  z),  and  if  the 
uncertainties  in  x,...,  z  are  independent  and  random,  then  the  uncertainty  in  the  computed  function  q  is 
given  by  (Taylor,  1982): 


(24) 
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We  treat  the  intensity  measurements  as  sadsfying  Equation  (24).  The  uncertainty  in  the  packing 
fraction,  Equadon  (22),  and  the  mass-thickness,  Equadon  (23),  result  from  uncertaindes  in  the  attenuation 
coefficients,  material  densides  and  the  x-ray  intensity'  measurements  which  are  independent  variables. 

The  largest  source  of  uncertainty  is  the  calculated  mass-attenuadon  coefficient  for  ammonia.  This 
uncertainty,  however,  is  not  random,  but  is  a  contribudon  to  the  systemadc  error.  The  contribudon  from 
the  ammoma  mass-attenuadon  coefficient  is  seven  dmes  more  significant  than  that  from  helium.  The 
accurate  knowledge  of  the  x-ray  energy  is  also  less  significant.  We  maintain  four  digits  of  accuracy  rather 
than  round  the  source  energy  up  to  60  keV  but  the  effect  of  lack  of  precision  in  the  x-ray  energy  is 
minimal.  For  example,  using  60.00  keV,  which  is  0.8%  higher  than  the  actual  x-ray  energy,  the 
denominator  in  Equadon  (23)  becomes  0.17400  resuldng  in  a  0.4%  decrease  in  thickness.  Ignoring  the 
systemadc  errors,  the  fracdonal  imcertainty  in  the  mass-thickness  is  a  fimcdon  of  intensides  only. 


5x.. 


_  Wy 
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where  the  uncertainty  in  intensides  has  been  evaluated  as  the  square  root  of  the  number  of  counts.  The 
last  step  is  to  show  explicidy  that  as  thickness  increases,  the  error  decreases  as  the  rado  goes  to 

infimty.  This  is  justified  because  idendcal  coundng  dmes  were  used.  The  quanddes  of  interest,  mass¬ 
thickness  and  volume  packing  fracdon,  describe  all  of  the  material  in  a  given  basket.  This  means  that  the 
intensides  used  for  our  calculadons  are  the  sums  of  the  21  intensides  from  each  scan  posidon.  This 
gready  improves  the  overall  accuracy'  as  point  to  point  fluctuadons  are  averaged  throughout  the  endre 
volume  of  a  cell.  Note  that  the  difference  between  the  average  of  the  logarithms  of  the  intensity  rados  and 
the  logarithm  of  the  average  of  the  intensity  rados  is  within  0.6%.  The  stadsdcal  error  on  the  final  values 
ranges  from  0.85%  to  1.3%.  The  systemadc  uncertainty  in  the  computer  tabulated  values  of  the 
attenuation  cross  secdons  are  not  well  established  experimentally,  but  theoredcally,  for  small  Z,  they  are 
exf)ected  to  be  in  the  2-5%  range  (McMaster,  V.,  et  al.^  1967). 

The  ammonia  samples  to  be  measured  were  placed  in  a  second  liquid  helium  cooled  cryostat  after 
being  transferred  from  the  actual  target  assembly  in  ESA.  The  posidon  of  a  target  cell  in  relation  to  the 
lab  table  outside  the  cryostat  was  determined  to  better  than  0.5  mm.  The  target  was  exposed  to  the  ''"Am 

source  for  five  rmnutes  per  scan  point.  The  configuradon  for  the  thickness  measurement  is  shown  in 
Figure  4  below. 
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Figure  4.  Schematic  of  target  thicknessdetermination  apparatus. 


One  important  difference  of  the  target  configuration  used  in  the  thickness  measurement  and  in  the 
electron  beam  runs  is  that  the  helium  bath  in  the  thickness  measurement  was  at  4  K  which  is  above  the 
lambda  point  of  helium  (2.18  K),  while  that  in  the  ESA  runs  was  below  the  lambda  point.  In  the  ESA 
electron  beam  runs  the  helium  density  is  approximately  13%  greater  than  during  the  thickness 
measurement.  Incorrectly  applying  the  greater  helium  density  of  0. 145  grams  per  cubic  cm  would 
increase  the  denominator  in  Equation  (22)  to  0.46847  for  NH3  and  0.46144  for  ND3,  resulting  in  a  1.8% 
increase  in  computed  packing  fraction.  The  packing  fraction  is  a  measurement  of  the  physical  dimensions 
of  the  ammonia,  irrespective  of  the  surrounding  helium.  The  packing  fraction  which  we  have  measured  is 
the  value  which  should  be  used. 

The  source  was  set  on  one  end  of  a  rectangular  steel  support  plate.  A  large  notch  was  cut  out  of  the 
side  of  the  support  plate  in  the  center  so  that  the  tailpiece  of  the  helium  dewar  could  be  set  in.  The 
support  plate  was  mounted  on  a  milling  machine  drive  that  allowed  a  Cartesian  coordinate  mapping  of  the 
target  face  by  adjusUng  the  support  plate  after  each  five  minute  scan  point  measurement.  At  several 
points  along  the  beam  path  there  are  lead  collimaUng  plates  with  holes  of  five  mm  diameter.  It  is  these 
holes  that  define  the  spots  that  overlap  and  cover  the  entire  target  face.  The  outermost  one  mm  ring  was 
not  scanned  by  the  electron  beam  in  ESA  and  so  was  omitted  from  the  thickness  mapping.  The  face  area 
of  a  target  window  is  covered,  with  minimal  overlap,  by  21  scan  points  as  shown  in  Figure  5  below. 
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Figure  5.  Pattern  of  21  scan  points  of  5  mm  diameter  each.  Target  face  is  25.4  mm  diameter. 

Because  of  the  presence  of  Cu/Ni  wire  coils  for  NMR  measurements  of  the  target  polarization,  the 
assignment  of  scan  point  centers  might  have  an  inordinate  influence  on  the  target  thickness  measurement. 
To  insure  that  this  was  not  the  case,  new  scan  point  locations  were  assigned  to  points  7,  9,  15  and  17. 
These  new  points  are  the  “primed”  points  T,  9’,  15’,  17’.  The  primed  scan  points  have  the  same 
horizontal  coordinates  as  the  unprimed  points  but  are  0.03  inches  closer  to  the  center  in  the  vertical 
direction.  Stick  number  one  has  no  prime  point  measurements  as  it  was  done  before  this  coverage  effect 
was  realized.  As  seen  in  Tables  6  through  12,  the  use  of  the  original  or  the  primed  points  results  in  very 
little  difference  in  the  final  thickness  results. 

The  detector  was  a  Canberra  Passivated  Implanted  Planar  Silicon  (PIPS)  with  associated  FET  input 
preamplifier  (Canberra  products  catalog,  1989).  Calibration  tests  were  performed  before  measurements 
began  to  ensure  that  the  detector  efficiency  and  resolution  were  constant  in  time  over  a  period  of  several 
hours.  The  output  signal  was  displayed  with  “Nucleus-IF’  personal  computer  analyzer  software 
(Oxford/Tennelec  Industries,  1990).  The  distance  from  the  Am  source  to  the  detector  was  28.5  cm. 
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IV.  DATA  AND  ANALYSIS 


Our  results  were  tabulated  using  “Excel”  spreadsheet  software  (Microsoft  Corp.,  1993).  There  are 
seven  data  sets,  each  consisting  of  a  top  cell  and  a  bottom  cell  entry.  All  the  ammonia  material  for  stick 
numbers  1  and  3  was  radiated  at  the  Stanford  High  Energy  Physics  Lab  60  MeV  linac,  but  the  ND3  for 
stick  number  2  was  from  Bates  Lab  while  the  NH3  was  radiated  at  the  Naval  Postgraduate  School  in 
Monterey.  The  data  was  divided  by  date.  The  data  for  stick  number  1  is  from  mid  January,  1994  and  did 
not  use  the  “prime”  raster  positions.  The  data  set  for  stick  number  2a  is  from  early  January  1994  and  did 
not  use  the  “prime”  raster  positions,  the  data  set  for  stick  number  2b  was  taken  later  that  month  and  also 
did  not  use  the  primed  scan  points.  The  data  for  stick  number  2b  prime  is  identical  to  number  2b  except 
for  using  the  values  at  the  primed  scan  point  locations.  The  data  set  for  stick  number  3a  is  from  early 
January,  1994,  while  the  data  set  for  stick  number  3b  was  taken  in  mid  February,  1994  after  final  exposure 
in  ESA.  The  data  for  stick  number  3b  prime  is  identical  to  number  3b  except  for  using  the  values  at  the 
primed  scan  point  locations.  The  data  analysis  software  (Oxford/Tennelec,  1990)  gives  count  intensities 
in  both  un-integrated  total  intensities  and  integrated  intensities  with  uniform  “background”  energies 
subtracted  out.  Our  calculations  use  the  first  value.  Our  assumption  in  including  the  “background”  from 
the  spectrum  was  that  it  was  directional  from  the  source  due  to  our  precautionary  shielding  of  the  detector 
and  collimating  of  the  beam.  Comparison  to  runs  of  much  longer  times  with  no  source  show  that  true 
atmospheric  background  was  completely  negligible  to  the  source  intensity  and  that  anything  getting  to  the 
detector  through  one  type  of  cell  was  getting  through  to  the  other  also.  The  intensity  calculations  of 
Equations  (19)  through  (23)  assume  total  counts  from  all  possible  sources. 


A.  ASSIGNMENT  OF  VALUES  TO  SCAN  POINTS 


The  columns  of  the  data  are  arranged  in  tables  as  follows;  “Scan  Point”  is  the  numbered  position  of  the 
five  mm  diameter  circle  of  the  pattern  of  2 1  used  to  cover  the  face  of  a  target.  The  next  two  columns  are 
grid  coordinates  of  the  relative  positions  of  the  scan  points.  The  fourth  column  is  the  integrated  intensity, 
including  “background,”  of  the  “empty”  target.  The  fifth  column  is  the  same  data  for  a  corresponding 
“full”  target,  again  including  the  background  since  we  are  subtracting  logarithms  and  so  must  include  the 
entire  sum  of  intensities  at  the  detector.  The  sixth  column  takes  the  logarithm  of  the  ratio  of  the  fourth  to 
the  fifth  column.  The  “Mass-Thik”  column  performs  Equation  (23)  by  dividing  the  sixth  column  by  the 
given  denominator  to  give  a  value  for  the  mass  thickness  in  grams  per  square  centimeter.  Column 
number  eight  gives  a  normalized  value  for  the  mass  thickness  by  dividing  each  point  by  the  average  of  the 
middle  nine  scan  points  (2  through  4,  7  through  9  and  15  through  17).  Column  number  nine  gives  the 
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packing  fraction  corresponding  to  the  thickness  value  in  the  seventh  column  by  performing  Equation  (22). 
The  last  column  gives  the  fractional  uncertainty  in  intensity  used  in  Equation  (25). 

The  bottom  two  entries  of  column  numbers  four  and  five  average  the  intensities  above  them  over  all 
21  points,  or  only  the  center  9,  while  for  the  last  five  column  calculations,  the  algebraic  operation  is 
performed  on  the  averages  derived  in  the  bottom  two  entries  of  column  number  six.  The  difference 
between  the  average  of  the  logarithms  and  the  logarithm  of  the  averages  is  negligible.  This  was 
determined  by  having  the  final  values  of  column  numbers  six  through  ten  calculated  by  averaging  over  the 
21  points  in  the  column  above  them  rather  than  operating  on  the  values  up  the  row  to  their  left.  The 
numbers  obtained  either  way  differ  by  less  than  0.2%. 

The  final  results  for  the  target  thickness  and  packing  fraction  measurement  are  displayed,  with  their 
random  errors,  in  Table  1  below. 


Stick  Dates  in  ESA  Sample  Ammonia  Thickness  Random  Packing  Random 


Number 

Beamline 

Material 

Cell  Type 

grams/cc 

Error 

Fraction 

Error 

1; 

1/13/94- 

1/23/94 

Stanford 

Stanford 

Top-NH3 

Bot-NH3 

1,729 

1.64 

0.009 

0.009 

0.629 

0.596 

0.003 

0.003 

2a; 

12/1/93- 

1/7/94 

Bates 

NFS 

Top-ND3 

Bot-NH3 

1.462 

1.747 

0.008 

0.009 

0.532 

0.551 

0.003 

0.003 

2b; 

12/1/93- 

1/7/94 

Bates 

NFS 

Top-ND3 

Bot-NH3 

1.53 

1.85 

0.008 

0.01 

0.556 

0.583 

0.003 

0.003 

2b'  ; 

12/1/93- 

1/7/94 

Bates 

NFS 

Top-ND3 

Bot-NH3 

1.5 

184 

0.008 

0.01 

0.546 

0.58 

0.003 

0.003 

3a; 

1/7/94- 

1/13/94 

HEFL 

HEFL 

Top-ND3 

Bot-NH3 

1.384 

1.643 

0.008 

0.009 

0.503 

0.519 

0.003 

0.003 

3b; 

1/23/94- 

completion 

HEFL 

HEFL 

Top-ND3 

Bot-NH3 

1.409 

1.701 

0.008 

0.009 

0.512 

0.537 

0.003 

0.003 

3b'; 

1/23/94- 

completion 

HEFL 

HEFL 

Top-ND3 

Bot-NH3 

1.41 

1.7 

0.008 

0.009 

0.512 

0.537 

0.003 

0.003 

Table  1.  Target  thickness  measurement  results 
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B.  THICKNESS  AS  A  FUNCTION  OF  HEIGHT 


An  important  question  about  the  experimental  results  is  whether  there  was  a  settling  effect  and  the 
average  thickness  decreases  with  height.  Table  2  shows  the  data  sets  arranged  by  vertical  elevation  on  the 
target  face,  averaged  over  horizontal  scan  points.  It  is  seen  that  the  effect  of  the  primed  scan  points  is 
minor  and  that  the  greatest  difference  due  to  using  them  over  the  unprimed  points  is  of  the  order  of  0.5% 
for  stick  number  3b  and  Stick  number  2b  ND3,  and  2.0%  for  the  NH3  in  stick  number  2b.  In  all  four 
cases,  the  correction  for  primed  scan  points  produced  lower  thickness  values.  Tables  3  through  12  and 
Figures  6  through  15  show  the  variation  of  target  thickness  as  a  function  of  vertical  position.  Note  that  in 
Figures  6  through  15,  the  error  bars  shown  are  the  ranee  of  the  horizontal  values  at  a  given  vertical 
height.  This  error  is  larger  for  the  central  positions  because  there  are  more  data  points  there  and  so  the 
range  of  values  spanned  by  the  average  is  greater.  All  of  the  graphs  appear  to  be  consistent  with  thickness 
decreasing  as  height  increases  except  stick  number  3b  prime  which  appears  relatively  flat  in  comparison. 

In  order  to  visualize  the  non-uniformity  of  the  target.  Figure  16  shows  the  grid  pattern  of  (unprimed) 
data  points  for  stick  number  1,  bottom  cell.  Figure  17  presents  a  three  dimensional  view  of  the  target 
thickness  for  the  bottom  cell  of  that  stick,  and  Figure  18  is  a  contour  plot  of  the  same  data.  The  primed 
scan  points  are  shown  in  Figure  19.  Figures  20  and  21  are  the  three  dimensional  and  contour  view  of  the 
bottom  cell  of  stick  number  3b  prime  and  show  a  more  uniform  thickness  which  explains  the  relative 
flatness  observed  in  Figure  14. 

The  measurement  of  the  target  thickness  was  simultaneously  performed  by  two  other  groups;  the  team 
of  Dr.  Ingo  Sick  and  Beni  Zihlmann  from  Basel  Switzerland  and  the  team  of  Dr.  Peter  Bosted  from  the 
American  University,  Washington,  D.C.  and  Todd  Averett  from  the  University  of  Virginia  at 
Charlottesville.  The  Basel  group  also  used  a  radiation  attenuation  measurement,  but  used  a  solid  carbon 
and  liquid  ammonia  target  produced  and  measured  in  Basel,  Switzerland  as  reference  targets  for  their 
ratios.  The  second  group  compared  event  rates  for  each  sample  after  the  fact  from  the  actual  experimental 
data  to  a  carbon  reference  target.  As  of  November,  1994,  there  was  still  some  error  between  the  three  sets 
of  data.  The  Averett/Bosted  data  is  approximately  in  the  middle  of  NPS  values,  which  are  all  lower,  and 
the  Basel  values,  which  are  all  higher.  The  sources  of  disagreement  and  systematic  error  have  not  been 
agreed  upon  as  yet  by  the  three  groups,  but  by  allowing  a  large  enough  uncertainty  to  encompass  the  range 
of  the  three  data  sets  for  ammonia  thickness,  the  final  uncertainty  in  the  Dilution  Factor  is  of  the  order  of 
(Bosted,  1994)  2%.  This  is  within  the  desired  accuracy  of  the  initial  proposal  (McCarthy,  J.,  et  ai,  1993). 
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APPENDIX.  GRAPHIC  REPRESENTATIONS  OF  DATA  SETS 


y-axis 

AvgThik 

;  PacFrac. 

y-axis 

AvgThik 

PacFrac. 

Stick  #1 

0.31 

1.508 

0.548164 

Stick  #1 

0.31 

1.474 

0.535805 

Top-NH3 

0.26 

1.072 

0.389676 

Bot-NH3 

0.26 

1.727 

0.627772 

Radiated 

0.17 

1.274 

0.463104 

Radiated 

0.17 

1.347 

0.48964 

at 

0 

1.762 

0.640494 

at 

0 

1.548 

0.562704 

HEPL, 

-0.17 

2.13 

0.774264 

HEPL, 

-0.17 

1.856 

0.674664 

Stanford 

-0.3 

2.26 

0.821519 

Stanford 

-0.3 

1.834 

0.666667 

-0.35 

2.505 

0.910578 

-0.35 

1.961 

0.712832 

Stick  #2a 

0.31 

1.046 

0.330177 

Stick  #2a 

0.31 

1.218 

0.442748 

Top-ND3 

0.26 

1.889 

0.596275 

Bot-NH3 

0.26 

0.6743 

0.245111 

Radiated 

0.17 

1.2 

0.378788 

Radiated 

0.17 

1.363 

0.495456 

at 

0 

1.841 

0.581124 

at 

0 

1.533 

0.557252 

Bates, 

-0.17 

2.448 

0.772727 

NPS, 

-0.17 

1.318 

0.479099 

MIT 

-0.3 

1.955 

0.617109 

Monterey 

-0.3 

1.534 

0.557615 

-0.35 

2.152 

0.679293 

-0.35 

2.648 

0.962559 

Stick  #2b 

0.31 

1.304 

0.411616 

Stick  #2b 

0.31 

0.7451 

0.270847 

Top-ND3 

0.26 

1.554 

0.49053 

Bot-NH3 

0.26 

1.15 

0.41803 

Radiated 

0.17 

1.659 

0.523674 

Radiated 

0.17 

1.068 

0.388222 

at 

0 

1.638 

0.517045 

at 

0 

1.761 

0.640131 

Bates, 

-0.17 

2.283 

0.720644 

NPS, 

-0.17 

1.818 

0.660851 

MIT 

-0.3 

2.018 

0.636995 

Monterey 

-0.3 

1.658 

0.60269 

-0.35 

2.463 

0.777462 

-0.35 

2.662 

0.967648 

Stick  #2b 

0.31 

1.304 

0.411616 

Stick  #2b 

0.31 

0.7451 

0.270847 

Top“ND3 

0.26 

1.554 

0.49053 

Bot-NH3 

0.26 

1.15 

0.41803 

Radiated 

0.14 

1.631 

0.514836 

Radiated 

0.14 

1.061 

0.385678 

at 

0 

1.638 

0.517045 

at 

0 

1.761 

0.640131 

Bates, 

-0.14 

2.264 

0.714646 

NPS, 

-0.14 

1.704 

0.619411 

MIT 

-0.3 

2.018 

0.636995 

Monterey 

-0.3 

1.658 

0.60269 

-0.35 

2.463 

0.777462 

-0.35 

2.662 

0.967648 

Stick  #3a 

0.31 

2.043 

0.644886 

Stick  #3a 

0.31 

1.45 

0.527081 

Top-ND3 

0.26 

1.891 

0.596907 

Bot-NH3 

0.26 

0.8156 

0.296474 

Radiated 

0.17 

1.161 

0.366477 

Radiated 

0.17 

1.901 

0.691021 

at 

0 

2.37 

0.748106 

at 

0 

1.276 

0.463831 

HEPL, 

-0.17 

1.184 

0.373737 

HEPL, 

-0.17 

1.055 

0.383497 

Stanford 

-0.3 

2.182 

0.688763 

Stanford 

-0.3 

1.821 

0.661941 

-0.35 

0.7641 

0.241193 

-0.35 

1.153 

0.41912 

Stick  #3a 

0.31 

0.9718 

0.306755 

Stick  #3b 

0.31 

0.9053 

0.32908 

Top-ND3 

0.26 

0.87 

0.274621 

Bot-NH3 

0.26 

1.138 

0.413668 

Radiated 

0.17 

1.52 

0.479798 

Radiated 

0.17 

1.393 

0.506361 

at 

0 

1.854 

0.585227 

at 

0 

1.462 

0.531443 

HEPL, 

-0.17 

1.877 

0.592487 

HEPL, 

-0.17 

1.484 

0.53944 

Stanford 

-0.3 

2.513 

0.793245 

Stanford 

-0.3 

1.521 

0.55289 

-0.35 

1.953 

0.616477 

-0.35 

1.696 

0.616503 

Stick  #3b 

0.31 

0.9718 

0.306755 

Stick  #3b 

0.31 

0.9053 

0.32908 

Top“ND3 

0.26 

0.87 

0.274621 

Bot-NH3 

0.26 

1.138 

0.413668 

Radiated 

0.14 

1.51 

0.476641 

Radiated 

0.14 

1.372 

0.498728 

at 

0 

1.854 

0.585227 

at 

0 

1.462 

0.531443 

HEPL, 

-0.14 

1.885 

0.595013 

HEPL, 

-0.14 

1.509 

0.548528 

Stanford 

-0.3 

2.513 

0.793245 

Stanford 

-0.3 

1.521 

0.55289 

-0.35 

1.953 

0.616477 

-0.35 

1.696 

0.616503 

Table  2.  Variation  of  data  sets  with  vertical  position. 
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Stickl 

-0.35 

2.552595 

2.457405 

2.505 

Top 

-0.3 

2.463802 

2.024597 

2.26 

15NH3 

-0.17 

2.471171 

1.51307 

2.13 

sample 

0 

2.562405 

1.20705 

1.762 

from 

0.17 

2.29905 

0.188916 

1.274 

HEPL 

0.26 

1.159399 

0.984245 

1.072 

0.31 

1.536652 

1.479348 

1.508 

Table  3.  Stick  #1  top  cell  variation  of  thickness  for  a  given  vertical  position. 


Stickl 

-0.35 

2.007672 

1.914328 

1,961 

Bottom 

1 

o 

2.017496 

1,85839 

1.834 

15NH3 

-0.17 

1.967024 

1.683543 

1.856 

Scimple 

0 

1.970723 

1.455216 

1.548 

from 

0.17 

2.017496 

0.924131 

1.347 

HEPL 

0.26 

1.837073 

0.909016 

1.727 

0.31 

1.508934 

1.439066 

1.474 

Table  4.  Stick  #1  bottom  cell  variation  of  thickness  for  a  given  vertical  position. 
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Stick2a 

-0.35 

2.202357 

2.101643 

2.152 

Top 

-0.3 

2.079346 

1.833116 

1.955 

15ND3 

-0.17 

2.82132 

2.123782 

2.448 

sample 

0 

2.815392 

0.872776 

1.841 

from 

0.17 

1.790826 

0.746651 

1.2 

Bates, 

0.26 

2.135189 

1.648595 

1.889 

MIT 

0.31 

1.078426 

1.013574 

1.046 

Table  5.  Stick  #2a  top  cell  variation  of  thickness  for  a  given  vertical  position. 


Figure  8.  Stick  #2a  top  cell  variation  of  thickness  for  a  given  vertical  position. 


Stick2a  -0.35 

2.698577 

2.597423 

2.648 

Bottom  -0-3 

1.688775 

1.381274 

1.534 

15NH3  -0.17 

1.673259 

1.022169 

1.318 

sample  0 

2.57016 

0.574774 

1.533 

from  0.17 

1.905008 

0.696056 

1.363 

NFS  0.26 

1.480258 

1.409742 

1.445 

0.31 

1.251008 

1.184992 

1.218 

Table  6.  Stick  #2a  bottom  cell  variation  of  thickness  for  a  given  vertical  position. 


Figure  9.  Stick  #2a  bottom  cell  variation  of  thickness  for  a  given  vertical  position. 


25 


Stick2b 

"0.35 

2.517679 

2.408321 

2.463 

Top 

-0.3 

2.074229 

1.961414 

2.017 

15ND3 

-0.14 

2.584666 

1.975662 

2.283 

primed 

0 

1.870592 

1.345529 

1.638 

sample 

0.14 

1.884318 

1.48413 

1.659 

from 

0.26 

1.655698 

1.454331 

1.544 

Bates 

0.31 

1.340512 

1.267488 

1.304 

Table  7.  Stick  #2b  top  cell  variation  of  thickness  for  a  given  vertical  position. 


Figure  10,  Stick  #2b  top  cell  variation  of  thickness  for  a  given  vertical  position. 


Stick2b 

-0.35 

2.712844 

2.611156 

2.662 

Bottom 

-0.3 

1.75257 

1.563566 

1.658 

15NH3 

-0.14 

2.383421 

1.467145 

1.818 

primed 

0 

2.115941 

1.387238 

1.761 

sample 

0.14 

1.774385 

0.400363 

1.068 

from 

o 

to 

2.191443 

0.141819 

1.15 

NFS 

0.31 

0.770359 

0.719841 

0.745 

Table  8.  Stick  #2b  bottom  cell  variation  of  thickness  for  a  given  vertical  position. 
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Stick3a 

-0*35 

0.792983 

0.735217 

0.764 

Top 

-0.3 

3.119259 

1.264312 

2.182 

15ND3 

-0.17 

1.723344 

0.6168 

1.184 

second 

0 

3.108664 

1.738526 

2.37 

sample 

0.17 

1.620907 

0.724594 

1.161 

from 

0.26 

2.066675 

1.716799 

1.891 

HEPL 

0.31 

2.088559 

1.961414 

2.043 

Table  9.  Stick  #3a  top  cell  variation  of  thickness  for  a  given  vertical  position. 


Figure  12.  Stick  #3a  top  cell  variation  of  thickness  for  a  given  vertical  position. 


Stick3a 

-0.35 

1.18563 

1.12037 

1.153 

Bottom 

-0.3 

2.259965 

1.391753 

1.821 

15NH3 

-0.17 

1.531927 

0.589776 

1.055 

second 

0 

1.713525 

1.054825 

1.276 

sample 

0.17 

2.215353 

1.053644 

1.901 

from 

0.26 

1.03618 

0.601625 

0.816 

HEPL 

0.31 

1.486105 

1.413895 

1.45 

Table  10.  Stick  #3a  bottom  cell  variation  of  thickness  for  a  given  vertical  position. 


Figure  13.  Stick  #3a  bottom  cell  variation  of  thickness  for  a  given  vertical  position. 
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Stick3b 

“0.35 

2.001434 

1.904566 

1.953 

Top 

“0.3 

2.722599 

2.304404 

2.513 

15ND3 

“0.14 

2.043904 

1.728763 

1.885 

primed 

0 

1.96026 

1.713856 

1.854 

second 

0.14 

1.737875 

1.320112 

1.51 

HEPL 

0.26 

1.010576 

0.732108 

0.87 

sample 

0.31 

1.001926 

0.941674 

0.972 

Table  11.  Stick  #3b  top  cell  variaUon  of  thickness  for  a  given  vertical  position. 


Figure  14.  Stick  #3b  top  cell  variation  of  thickness  for  a  given  vertical  position. 


Stick3b 

-0.35 

1.736704 

1.655296 

1.696 

Bottom 

-0.3 

1.678579 

1.366115 

1.521 

15NH3 

-0.14 

1.68501 

1.396057 

1.509 

primed 

0 

1.679196 

1.210389 

1.462 

second 

0.14 

1.538471 

1.220547 

1.372 

HEPL 

0.26 

1.240313 

1.038756 

1.138 

sample 

0.31 

0.932784 

0.877217 

0,905 

Table  12.  Stick  #3b  bottom  cell  variaUon  of  thickness  for  a  given  verUcal  posiUon. 


Figure  15.  SUck  #3b  bottom  cell  variation  of  thickness  for  a  given  vertical  position 
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1/100  of  on  inch  vertical  extent 


Figure  16.  Centers  of  the  “unprimed”  scan  points  #’s  1-21. 
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Figure  17.  Thickness  variaUon  for  suck  #1,  bottom,  Bates  NHj  target.  Data  presented  in  Table  1. 
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10 


20 

1/100  of 


30  40  50 

on  inch  horizontal  extent 


Figure  18.  Thickness  contour  for  target  shown  in  Figure  6. 


Figure  20.  Thickness  variation  for  stick  #3b  prime,  bottom,  HEPL  NH3  target.  Data  presented  in  Table  9. 
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l/lOO  of  on  inch  vertical  extent 


Figure  21.  Thickness  contour  for  target  shown  in  Figure  9. 
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